SARS-CoV-2 intra-host diversity, antibody response, and disease severity after reinfection by the variant of concern Gamma in Brazil.

The rapid spread of the SARS-CoV-2 Variant of Concern (VOC) Gamma in Amazonas during early 2021 fueled a second large COVID-19 epidemic wave and raised concern about the potential role of reinfections. Very few cases of reinfection associated with the VOC Gamma have been reported to date, and their potential impact on clinical, immunological, and virological parameters remains largely unexplored. Here we describe 25 cases of SARS-CoV-2 reinfection in Brazil. SARS-CoV-2 genomic analysis confirmed that individuals were primo-infected with distinct viral lineages between March and December 2020 (B.1.1, B.1.1.28, B.1.1.33, B.1.195, and P.2) and reinfected with the VOC Gamma between 3 to 12 months after primo-infection. We found a similar mean cycle threshold (Ct) value and limited intra-host viral diversity in both primo-infection and reinfection samples. Sera of 14 patients tested 10-75 days after reinfection displayed detectable neutralizing antibodies (NAb) titers against SARS-CoV-2 variants that circulated before (B.1.*), during (Gamma), and after (Delta and Omicron) the second epidemic wave in Brazil. All individuals had milder or no symptoms after reinfection, and none required hospitalization. These findings demonstrate that individuals reinfected with the VOC Gamma may display relatively high RNA viral loads at the upper respiratory tract after reinfection, thus contributing to onward viral transmissions. Despite this, our study points to a low overall risk of severe Gamma reinfections, supporting that the abrupt increase in hospital admissions and deaths observed in Amazonas and other Brazilian states during the Gamma wave was mostly driven by primary infections. Our findings also indicate that most individuals analyzed developed a high anti-SARS-CoV-2 NAb response after reinfection that may provide some protection against reinfection or disease by different SARS-CoV-2 variants.

Comparative epidemic expansion of SARS-CoV-2 variants Delta and Omicron in the Brazilian State of Amazonas.

The SARS-CoV-2 variants of concern (VOCs) Delta and Omicron spread globally during mid and late 2021, respectively. In this study, we compare the dissemination dynamics of these VOCs in the Amazonas state, one of Brazil's most heavily affected regions. We sequenced the virus genome from 4128 patients collected in Amazonas between July 1st, 2021, and January 31st, 2022, and investigated the viral dynamics using a phylodynamic approach. The VOCs Delta and Omicron BA.1 displayed similar patterns of phylogeographic spread but different epidemic dynamics. The replacement of Gamma by Delta was gradual and occurred without an upsurge of COVID-19 cases, while the rise of Omicron BA.1 was extremely fast and fueled a sharp increase in cases. Thus, the dissemination dynamics and population-level impact of new SARS-CoV-2 variants introduced in the Amazonian population after mid-2021, a setting with high levels of acquired immunity, greatly vary according to their viral phenotype.

Molecular characterization of a new SARS-CoV-2 recombinant cluster XAG identified in Brazil.

Recombination events have been described in the Coronaviridae family. Since the beginning of the SARS-CoV-2 pandemic, a variable degree of selection pressure has acted upon the virus, generating new strains with increased fitness in terms of viral transmission and antibody scape. Most of the SC2 variants of concern (VOC) detected so far carry a combination of key amino acid changes and indels. Recombination may also reshuffle existing genetic profiles of distinct strains, potentially giving origin to recombinant strains with altered phenotypes. However, co-infection and recombination events are challenging to detect and require in-depth curation of assembled genomes and sequencing reds. Here, we present the molecular characterization of a new SARS-CoV-2 recombinant between BA.1.1 and BA.2.23 Omicron lineages identified in Brazil. We characterized four mutations that had not been previously described in any of the recombinants already identified worldwide and described the likely breaking points. Moreover, through phylogenetic analysis, we showed that the newly named XAG lineage groups in a highly supported monophyletic clade confirmed its common evolutionary history from parental Omicron lineages and other recombinants already described. These observations were only possible thanks to the joint effort of bioinformatics tools auxiliary in genomic surveillance and the manual curation of experienced personnel, demonstrating the importance of genetic, and bioinformatic knowledge in genomics.

Real-life evaluation of a rapid antigen test (DPP® SARS-CoV-2 antigen) for COVID-19 diagnosis of primary healthcare patients, in the context of the omicron-dominant wave in Brazil.

OBJECTIVES: We aimed to investigate the real-life performance of the rapid antigen test (RAT) in the context of a primary healthcare setting, including symptomatic and asymptomatic individuals that sought diagnostic during a Omicron infection wave. METHODS: We prospectively accessed the performance of the DPP® SARS-CoV-2 Antigen test in the context of an omicron-dominant real-life setting. We evaluated 347 unselected individuals (all-comers) from a public testing center in Brazil, performing the RAT diagnosis at point-of-care with fresh samples. The combinatory result from two distinct RT-qPCR methods was employed as reference and 13 samples with discordant PCR results were excluded. RESULTS: The assessment of the rapid test in 67 PCR-positive and 265 negative samples revealed an overall sensitivity of 80.5% (CI95% = 69.1 - 89.2%), specificity of 99.2% (CI95% = 97.3 - 99.1%) and positive/negative predictive values higher than 95%. However, we observed that the sensitivity was dependent on the viral load (sensitivity in Ct<31 = 93.7%, CI = 82.8 - 98.7%; Ct>31 = 47.4%, CI = 24.4 - 71.1%). The positive samples evaluated in the study were Omicron (BA.1/BA.1.1) by whole-genome sequencing (n=40) and multiplex RT-qPCR (n=17). CONCLUSIONS: Altogether, the data obtained from a real-life prospective cohort supports that the RAT sensitivity for Omicron remains high and underscores the reliability of the test for COVID-19 diagnosis in settings with high disease prevalence and limited PCR testing capability.

Molecular characterization of a new SARS-CoV-2 recombinant cluster XAG identified in Brazil

Recombination events have been described in the Coronaviridae family. Since the beginning of the SARS-CoV-2 pandemic, a variable degree of selection pressure has acted upon the virus, generating new strains with increased fitness in terms of viral transmission and antibody scape. Most of the SC2 variants of concern (VOC) detected so far carry a combination of key amino acid changes and indels. Recombination may also reshuffle existing genetic profiles of distinct strains, potentially giving origin to recombinant strains with altered phenotypes. However, co-infection and recombination events are challenging to detect and require in-depth curation of assembled genomes and sequencing reds. Here, we present the molecular characterization of a new SARS-CoV-2 recombinant between BA.1.1 and BA.2.23 Omicron lineages identified in Brazil. We characterized four mutations that had not been previously described in any of the recombinants already identified worldwide and described the likely breaking points. Moreover, through phylogenetic analysis, we showed that the newly named XAG lineage groups in a highly supported monophyletic clade confirmed its common evolutionary history from parental Omicron lineages and other recombinants already described. These observations were only possible thanks to the joint effort of bioinformatics tools auxiliary in genomic surveillance and the manual curation of experienced personnel, demonstrating the importance of genetic, and bioinformatic knowledge in genomics.

Phylogenetic-based inference reveals distinct transmission dynamics of SARS-CoV-2 lineages Gamma and P.2 in Brazil.

The COVID-19 epidemic in Brazil experienced two major lineage replacements until mid-2021. The first was driven by lineage P.2, in late 2020, and the second by lineage Gamma, in early 2021. To understand how these SARS-CoV-2 lineages spread in Brazil, we analyzed 11,724 genomes collected throughout the country between September 2020 and April 2021. Our findings indicate that lineage P.2 probably emerged in July 2020 in the Rio de Janeiro state and Gamma in November 2020 in the Amazonas state. Both states were the main hubs of viral disseminations to other Brazilian locations. We estimate that Gamma was 1.56-3.06 times more transmissible than P.2 in Rio de Janeiro and that the median effective reproductive number (Re) of Gamma varied according to the geographic context (Re = 1.59-3.55). In summary, our findings support that lineage Gamma was more transmissible and spread faster than P.2 in Brazil.

Unusual SARS-CoV-2 intrahost diversity reveals lineage superinfection.

Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has infected almost 200 million people worldwide by July 2021 and the pandemic has been characterized by infection waves of viral lineages showing distinct fitness profiles. The simultaneous infection of a single individual by two distinct SARS-CoV-2 lineages may impact COVID-19 disease progression and provides a window of opportunity for viral recombination and the emergence of new lineages with differential phenotype. Several hundred SARS-CoV-2 lineages are currently well phylogenetically defined, but two main factors have precluded major coinfection/codetection and recombination analysis thus far: (i) the low diversity of SARS-CoV-2 lineages during the first year of the pandemic, which limited the identification of lineage defining mutations necessary to distinguish coinfecting/recombining viral lineages; and the (ii) limited availability of raw sequencing data where abundance and distribution of intrasample/intrahost variability can be accessed. Here, we assembled a large sequencing dataset from Brazilian samples covering a period of 18 May 2020 to 30 April 2021 and probed it for unexpected patterns of high intrasample/intrahost variability. This approach enabled us to detect nine cases of SARS-CoV-2 coinfection with well characterized lineage-defining mutations, representing 0.61 % of all samples investigated. In addition, we matched these SARS-CoV-2 coinfections with spatio-temporal epidemiological data confirming its plausibility with the cocirculating lineages at the timeframe investigated. Our data suggests that coinfection with distinct SARS-CoV-2 lineages is a rare phenomenon, although it is certainly a lower bound estimate considering the difficulty to detect coinfections with very similar SARS-CoV-2 lineages and the low number of samples sequenced from the total number of infections.

Identification of a novel SARS-CoV-2 P.1 sub-lineage in Brazil provides new insights about the mechanisms of emergence of variants of concern.

One of the most remarkable severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC) features is the significant number of mutations they acquired. However, the specific factors that drove the emergence of such variants since the second half of 2020 are not fully resolved. In this study, we describe a new SARS-CoV-2 P.1 sub-lineage circulating in Brazil, denoted here as Gamma-like-II, that as well as the previously described lineage Gamma-like-I shares several lineage-defining mutations with the VOC Gamma. Reconstructions of ancestor sequences support that most lineage-defining mutations of the Spike (S) protein, including those at the receptor-binding domain (RBD), accumulated at the first P.1 ancestor. In contrast, mutations outside the S protein were mostly fixed at subsequent steps. Our evolutionary analyses estimate that P.1-ancestral strains carrying RBD mutations of concern probably circulated cryptically in the Amazonas for several months before the emergence of the VOC Gamma. Unlike the VOC Gamma, the other P.1 sub-lineages displayed a much more restricted dissemination and accounted for a low fraction (<2 per cent) of SARS-CoV-2 infections in Brazil in 2021. The stepwise diversification of lineage P.1 through multiple inter-host transmissions is consistent with the hypothesis that partial immunity acquired from natural SARS-CoV-2 infections in heavily affected regions might have been a major driving force behind the natural selection of some VOCs. The lag time between the emergence of the P.1 ancestor and the expansion of the VOC Gamma and the divergent epidemic trajectories of P.1 sub-lineages support a complex interplay between the emergence of mutations of concern and viral spread in Brazil.

Spread of Gamma (P.1) Sub-Lineages Carrying Spike Mutations Close to the Furin Cleavage Site and Deletions in the N-Terminal Domain Drives Ongoing Transmission of SARS-CoV-2 in Amazonas, Brazil.

The Amazonas was one of the most heavily affected Brazilian states by the COVID-19 epidemic. Despite a large number of infected people, particularly during the second wave associated with the spread of the Variant of Concern (VOC) Gamma (lineage P.1), SARS-CoV-2 continues to circulate in the Amazonas. To understand how SARS-CoV-2 persisted in a human population with a high immunity barrier, we generated 1,188 SARS-CoV-2 whole-genome sequences from individuals diagnosed in the Amazonas state from 1st January to 6th July 2021, of which 38 were vaccine breakthrough infections. Our study reveals a sharp increase in the relative prevalence of Gamma plus (P.1+) variants, designated Pango Lineages P.1.3 to P.1.6, harboring two types of additional Spike changes: deletions in the N-terminal (NTD) domain (particularly Δ144 or Δ141-144) associated with resistance to anti-NTD neutralizing antibodies or mutations at the S1/S2 junction (N679K or P681H) that probably enhance the binding affinity to the furin cleavage site, as suggested by our molecular dynamics simulations. As lineages P.1.4 (S:N679K) and P.1.6 (S:P681H) expanded (Re > 1) from March to July 2021, the lineage P.1 declined (Re < 1) and the median Ct value of SARS-CoV-2 positive cases in Amazonas significantly decreases. Still, we did not find an increased incidence of P.1+ variants among breakthrough cases of fully vaccinated patients (71%) in comparison to unvaccinated individuals (93%). This evidence supports that the ongoing endemic transmission of SARS-CoV-2 in the Amazonas is driven by the spread of new local Gamma/P.1 sublineages that are more transmissible, although not more efficient to evade vaccine-elicited immunity than the parental VOC. Finally, as SARS-CoV-2 continues to spread in human populations with a declining density of susceptible hosts, the risk of selecting more infectious variants or antibody evasion mutations is expected to increase. IMPORTANCE The continuous evolution of SARS-CoV-2 is an expected phenomenon that will continue to happen due to the high number of cases worldwide. The present study analyzed how a Variant of Concern (VOC) could still circulate in a population hardly affected by two COVID-19 waves and with vaccination in progress. Our results showed that the answer behind that was a new generation of Gamma-like viruses, which emerged locally carrying mutations that made it more transmissible and more capable of spreading, partially evading prior immunity triggered by natural infections or vaccines. With thousands of new cases daily, the current pandemics scenario suggests that SARS-CoV-2 will continue to evolve and efforts to reduce the number of infected subjects, including global equitable access to COVID-19 vaccines, are mandatory. Thus, until the end of pandemics, the SARS-CoV-2 genomic surveillance will be an essential tool to better understand the drivers of the viral evolutionary process.

ViralFlow: A Versatile Automated Workflow for SARS-CoV-2 Genome Assembly, Lineage Assignment, Mutations and Intrahost Variant Detection.

The COVID-19 pandemic is driven by Severe Acute Respiratory Syndrome coronavirus 2 (SARS-CoV-2) that emerged in 2019 and quickly spread worldwide. Genomic surveillance has become the gold standard methodology used to monitor and study this fast-spreading virus and its constantly emerging lineages. The current deluge of SARS-CoV-2 genomic data generated worldwide has put additional pressure on the urgent need for streamlined bioinformatics workflows. Here, we describe a workflow developed by our group to process and analyze large-scale SARS-CoV-2 Illumina amplicon sequencing data. This workflow automates all steps of SARS-CoV-2 reference-based genomic analysis: data processing, genome assembly, PANGO lineage assignment, mutation analysis and the screening of intrahost variants. The pipeline is capable of processing a batch of around 100 samples in less than half an hour on a personal laptop or in less than five minutes on a server with 50 threads. The workflow presented here is available through Docker or Singularity images, allowing for implementation on laptops for small-scale analyses or on high processing capacity servers or clusters. Moreover, the low requirements for memory and CPU cores and the standardized results provided by ViralFlow highlight it as a versatile tool for SARS-CoV-2 genomic analysis.

Identification of a novel SARS-CoV-2 P.1 sub-lineage in Brazil provides new insights about the mechanisms of emergence of variants of concern

One of the most remarkable severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern (VOC) features is the significant number of mutations they acquired. However, the specific factors that drove the emergence of such variants since the second half of 2020 are not fully resolved. In this study, we describe a new SARS-CoV-2 P.1 sub-lineage circulating in Brazil, denoted here as Gamma-like-II, that as well as the previously described lineage Gamma-like-I shares several lineage-defining mutations with the VOC Gamma. Reconstructions of ancestor sequences support that most lineage-defining mutations of the Spike (S) protein, including those at the receptor-binding domain (RBD), accumulated at the first P.1 ancestor. In contrast, mutations outside the S protein were mostly fixed at subsequent steps. Our evolutionary analyses estimate that P.1-ancestral strains carrying RBD mutations of concern probably circulated cryptically in the Amazonas for several months before the emergence of the VOC Gamma. Unlike the VOC Gamma, the other P.1 sub-lineages displayed a much more restricted dissemination and accounted for a low fraction (<2 per cent) of SARS-CoV-2 infections in Brazil in 2021. The stepwise diversification of lineage P.1 through multiple inter-host transmissions is consistent with the hypothesis that partial immunity acquired from natural SARS-CoV-2 infections in heavily affected regions might have been a major driving force behind the natural selection of some VOCs. The lag time between the emergence of the P.1 ancestor and the expansion of the VOC Gamma and the divergent epidemic trajectories of P.1 sub-lineages support a complex interplay between the emergence of mutations of concern and viral spread in Brazil.

Optimization and Clinical Validation of Colorimetric Reverse Transcription Loop-Mediated Isothermal Amplification, a Fast, Highly Sensitive and Specific COVID-19 Molecular Diagnostic Tool That Is Robust to Detect SARS-CoV-2 Variants of Concern.

The coronavirus disease 2019 (COVID-19) pandemic unfolded due to the widespread severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) transmission reinforced the urgent need for affordable molecular diagnostic alternative methods for massive testing screening. We present the clinical validation of a pH-dependent colorimetric reverse transcription loop-mediated isothermal amplification (RT-LAMP) for SARS-CoV-2 detection. The method revealed a limit of detection of 19.3 ± 2.7 viral genomic copies/µL when using RNA extracted samples obtained from nasopharyngeal swabs collected in guanidine-containing viral transport medium. Typical RT-LAMP reactions were performed at 65°C for 30 min. When compared to reverse transcriptase-quantitative polymerase chain reaction (RT-qPCR), up to cycle-threshold (Ct) value 32, RT-LAMP presented 98% [95% confidence interval (CI) = 95.3-99.5%] sensitivity and 100% (95% CI = 94.5-100%) specificity for SARS-CoV-2 RNA detection targeting E and N genes. No cross-reactivity was detected when testing other non-SARS-CoV virus, confirming high specificity. The test is compatible with primary RNA extraction-free samples. We also demonstrated that colorimetric RT-LAMP can detect SARS-CoV-2 variants of concern and variants of interest, such as variants occurring in Brazil named gamma (P.1), zeta (P.2), delta (B.1.617.2), B.1.1.374, and B.1.1.371. The method meets point-of-care requirements and can be deployed in the field for high-throughput COVID-19 testing campaigns, especially in countries where COVID-19 testing efforts are far from ideal to tackle the pandemics. Although RT-qPCR is considered the gold standard for SARS-CoV-2 RNA detection, it requires expensive equipment, infrastructure, and highly trained personnel. In contrast, RT-LAMP emerges as an affordable, inexpensive, and simple alternative for SARS-CoV-2 molecular detection that can be applied to massive COVID-19 testing campaigns and save lives.

COVID-19 in Amazonas, Brazil, was driven by the persistence of endemic lineages and P.1 emergence.

The northern state of Amazonas is among the regions in Brazil most heavily affected by the COVID-19 epidemic and has experienced two exponentially growing waves, in early and late 2020. Through a genomic epidemiology study based on 250 severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genomes from different Amazonas municipalities sampled between March 2020 and January 2021, we reveal that the first exponential growth phase was driven mostly by the dissemination of lineage B.1.195, which was gradually replaced by lineage B.1.1.28 between May and June 2020. The second wave coincides with the emergence of the variant of concern (VOC) P.1, which evolved from a local B.1.1.28 clade in late November 2020 and replaced the parental lineage in <2 months. Our findings support the conclusion that successive lineage replacements in Amazonas were driven by a complex combination of variable levels of social distancing measures and the emergence of a more transmissible VOC P.1 virus. These data provide insights to understanding the mechanisms underlying the COVID-19 epidemic waves and the risk of dissemination of SARS-CoV-2 VOC P.1 in Brazil and, potentially, worldwide.

A Sanger-based approach for scaling up screening of SARS-CoV-2 variants of interest and concern.

The global spread of new SARS-CoV-2 variants of concern underscore an urgent need of simple deployed molecular tools that can differentiate these lineages. Several tools and protocols have been shared since the beginning of the COVID-19 pandemic, but they need to be timely adapted to cope with SARS-CoV-2 evolution. Although whole-genome sequencing (WGS) of the virus genetic material has been widely used, it still presents practical difficulties such as high cost, shortage of available reagents in the global market, need of a specialized laboratorial infrastructure and well-trained staff. These limitations result in SARS-CoV-2 surveillance blackouts across several countries. Here we propose a rapid and accessible protocol based on Sanger sequencing of a single PCR fragment that is able to identify and discriminate all SARS-CoV-2 variants of concern (VOCs) identified so far, according to each characteristic mutational profile at the Spike-RBD region (K417N/T, E484K, N501Y, A570D). Twelve COVID-19 samples from Brazilian patients were evaluated for both WGS and Sanger sequencing: three P.2, two P.1, six B.1.1 and one B.1.1.117 lineage. All results from the Sanger sequencing method perfectly matched the mutational profile of VOCs and non-VOCs RBD's characterized by WGS. In summary, this approach allows a much broader network of laboratories to perform molecular surveillance of SARS-CoV-2 VOCs and report results within a shorter time frame, which is of utmost importance in the context of rapid public health decisions in a fast evolving worldwide pandemic.

Evolutionary Dynamics and Dissemination Pattern of the SARS-CoV-2 Lineage B.1.1.33 During the Early Pandemic Phase in Brazil.

A previous study demonstrates that most of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) Brazilian strains fell in three local clades that were introduced from Europe around late February 2020. Here we investigated in more detail the origin of the major and most widely disseminated SARS-CoV-2 Brazilian lineage B.1.1.33. We recovered 190 whole viral genomes collected from 13 Brazilian states from February 29 to April 31, 2020 and combined them with other B.1.1 genomes collected globally. Our genomic survey confirms that lineage B.1.1.33 is responsible for a variable fraction of the community viral transmissions in Brazilian states, ranging from 2% of all SARS-CoV-2 genomes from Pernambuco to 80% of those from Rio de Janeiro. We detected a moderate prevalence (5-18%) of lineage B.1.1.33 in some South American countries and a very low prevalence (<1%) in North America, Europe, and Oceania. Our study reveals that lineage B.1.1.33 evolved from an ancestral clade, here designated B.1.1.33-like, that carries one of the two B.1.1.33 synapomorphic mutations. The B.1.1.33-like lineage may have been introduced from Europe or arose in Brazil in early February 2020 and a few weeks later gave origin to the lineage B.1.1.33. These SARS-CoV-2 lineages probably circulated during February 2020 and reached all Brazilian regions and multiple countries around the world by mid-March, before the implementation of air travel restrictions in Brazil. Our phylodynamic analysis also indicates that public health interventions were partially effective to control the expansion of lineage B.1.1.33 in Rio de Janeiro because its median effective reproductive number (R e ) was drastically reduced by about 66% during March 2020, but failed to bring it to below one. Continuous genomic surveillance of lineage B.1.1.33 might provide valuable information about epidemic dynamics and the effectiveness of public health interventions in some Brazilian states.

Multiple Introductions Followed by Ongoing Community Spread of SARS-CoV-2 at One of the Largest Metropolitan Areas of Northeast Brazil.

Multiple epicenters of the SARS-CoV-2 pandemic have emerged since the first pneumonia cases in Wuhan, China, such as Italy, USA, and Brazil. Brazil is the third-most affected country worldwide, but genomic sequences of SARS-CoV-2 strains are mostly restricted to states from the Southeast region. Pernambuco state, located in the Northeast region, is the sixth most affected Brazilian state, but very few genomic sequences from the strains circulating in this region are available. We sequenced 101 strains of SARS-CoV-2 from patients presenting Covid-19 symptoms that reside in Pernambuco. Phylogenetic reconstructions revealed that all genomes belong to the B lineage and most of the samples (88%) were classified as lineage B.1.1. We detected multiple viral introductions from abroad (likely from Europe) as well as six local B.1.1 clades composed by Pernambuco only strains. Local clades comprise sequences from the capital city (Recife) and other country-side cities, corroborating the community spread between different municipalities of the state. These findings demonstrate that different from Southeastern Brazilian states where the epidemics were majorly driven by one dominant lineage (B.1.1.28 or B.1.1.33), the early epidemic phase at the Pernambuco state was driven by multiple B.1.1 lineages seeded through both national and international traveling.